Determination of Penetration Depth of Transverse Spin Current in Ferromagnetic 

Metals by Spin Pumping 
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Spin pumping in nonmagnetic/ferromagnetic metal multilayers is studied both theoretically and 
experimentally. We show that the line widths of the ferromagnetic resonance (FMR) spectrum 
depend on the thickness of the ferromagnetic metal layers, which must not be in resonance with the 
oscillating magnetic field. We also show that the penetration depths of the transverse spin current 
in ferromagnetic metals can be determined by analyzing the line widths of their FMR spectra. 
The obtained penetration depths in NiFe, CoFe and CoFeB were 3.7 [nm], 2.5 [nm] and 12.0 [nm], 
respectively. 

PACS numbers; 72.25.Mk, 75.70.Cn, 76.50.+g, 76.60.Es 
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The field of current-driven magnetization dynamics 
(CDMD) has drawn enormous attention because of its 
potential applications to non-volatile magnetic random 
access memory and microwave devices. CDMD is also 
important from a scientific point of view since it pro- 
vides much information about non-equilibrium dynamics 
of the magnetization and physics of spin transport and 
spin relaxation. The concept of CDMD was first pro- 
posed by Slonczewski [l| and independently by Berger 
[2] in 1996. In the last decade much effort has been de- 
voted to studying the physics and applications of CDMD 
both theoretically and experimentally [1, |3| • 

One of the most important quantities in CDMD is the 
penetration depth of the transverse spin current At , over 
which spin transfer torque is exerted for the magnetiza- 
tion of the free layer. However, there is a controversial 
issue regarding the penetration depth of the transverse 
spin current. One argument is based on the ballistic the- 
ory of electron transport, and its At = 7r/|fcp — /cp|, which 
is on the order of the lattice constant in conventional fer- 
romagnets such as Fe, Co, Ni, and their alloys [1, Q. 
The other argument is based on the Boltzmann theory 
of electron transport, and its At is on the order of a few 
nm 0] ■ Urazhdin et al. analyzed the CPP-GMR of non- 
collinear magnetic multilayers using the extended two- 
series-resistance model and concluded that At=0.8 [nm] 
for permalloy 8]. On the other hand, Chen et al. ana- 
lyzed the critical current of the CDMD in the Co/Cu/Co 
trilayer system and concluded that At =3.0 [nm] for Co 
i- 

The inverse process of CDMD is spin pumping, where 
spin current is generated by precession of magnetiza- 
tion in the ferromagnetic layer p^ . Enhancement of 
the Gilbert damping constant due to spin pumping has 
been extensively studied, and spin diffusion lengths, i.e., 
penetration depths of spin current in nonmagnetic met- 
als, have been obtained by analyzing the dependence of 
the enhancement of the Gilbert damping constant on 
the thickness of the nonmagnetic metal layer. In spin 
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FIG. 1: Schematic illustration of a nonmagnetic/ferromag- 
netic metal five-layer system. The magnetization of the Fi 
layer (mi) precesses around the z-sxis with angle 9. The 
magnetization of the F2 layer (m2) is fixed with the z-axis. 
The precession of the magnetization in the Fi layer pumps 
the spin current 1^"™^. The pumped spin current creates 
spin accumulation in the other layers, and the spin accumula- 
tion induces a backflow of spin current I^^^ across each N/F 
interface. 



pumping, the direction of the magnetization vector of 
the pumped spin current is perpendicular to the direc- 
tion of the precessing magnetization vector [To| . Let us 
consider the nonmagnetic/ferromagnetic metal five-layer 
system shown in Fig. [TJ Since the magnetization vec- 
tor of the pumped spin current IP"™p is perpendicular 
to the magnetization vector mi of the Fi layer and the 
precession angle 9 is very small (about 1 [deg]) in con- 
ventional FMR experiments, the dominant component of 
the pumped spin current is perpendicular to the magne- 
tization vector m2 of the F2 layer. Therefore, it would be 
possible to determine the penetration depth of the trans- 
verse spin current in the F2 layer if we could analyze the 
dependence of the enhancement of the Gilbert damping 
constant on the thickness of the F2 layer. 

In this letter, we study spin pumping in N1/F1/N2/ 
F2/N3 five-layer systems shown in Fig. [T] both theo- 
retically and experimentally. We extend Tserkovnyak's 
theory of spin pumping by taking into account the fi- 
nite penetration depth of the transverse spin current and 
show that the enhancement of the Gilbert damping con- 
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stant due to spin pumping depends on the ratio of the 
penetration depth At and the thickness d2 of the F2 layer. 
The motion of the magnetization vector in the F2 layer 
is not in resonance with an oscillating magnetic field; 
hence, the F2 layer plays the role of spin absorber. We 
also perform FMR experiments in Cu/CoFe/Cu/Py/Cu, 
Cu/Py/Cu/CoFe/Cu and Cu/CoFe/Cu/CoFeB/Cu five- 
layer systems and measure line widths AB of energy ab- 
sorption spectra, which are closely related to the Gilbert 
damping constants. The abbreviations CoFe, CoFeB, 
and Py hereafter refer to Co75Fe25, (Co5oFe5o)8oB2o and 
Ni8oFe2o, respectively. Analyzing the dependence of the 
line width on the thickness of the Py, CoFe and CoFeB 
layers that are not in resonance, we showed that the pen- 
etration depths of the transverse spin current in the Py, 
CoFe and CoFeB layers are 3.7 [nm], 2.5 [nm] and 12.0 
[nm], respectively. 

Let us begin with an introduction to the theory of spin 
pumping in the nonmagnetic/ferromagnetic metal five- 
layer system shown in Fig. [1] with a finite penetration 
depth of the transverse spin current. The pumped spin 
current generated by precession of the magnetization mi 
of the Fi layer is given by 
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where h is the Dirac constant and g^^l^^ is the real (imag- 
inary) part of the mixing conductance [13. The pumped 
spin current creates spin accumulation in the other layers, 
and the spin accumulation induces a backflow of spin cur- 
rent across each N/F interface. Although the backflow 
is obtained from circuit theory [l3|, the penetration 
depth of the transverse spin current At is assumed to be 
zero in this theory. Since we are interested in the effect 
of the penetration depth of the transverse spin current 
on spin pumping, we explicitly consider the diffusion pro- 
cess of transverse spin accumulation in the ferromagnetic 
layer. The backflow of spin current flowing from the 
layer to the Ffc layer is expressed as 
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where (;TT(ii) jg i^jig gpjji (down) conductance, 
is the real (imaginary) part of the transmission mixing 
conductance at the Ffc/N^ interface and /iNi is the spin 
accumulation in the layer The longitudinal spin 
accumulation in the Ffc layer is denoted by Mp^. The 
last two terms express contributions from transverse spin 
accumulation fij.^ in the Ffc layer. 

The spin accumulation in a ferromagnetic layer is de- 
fined by the non-equilibrium distribution matrix at a 
given energy e, /(e) = /oi + / • <t 0, where / = 



/i:*! -|-/j,t2 -|-/zm. Here we introduce the orthogonal unit 
vectors in spin space (ti,t2,m). The non-equilibrium 
charge distribution is represented by /o = (/^ + /^)/2. 
On the other hand, fz = (/^ — f^)/'2, is the difference in 
non-equilibrium distribution between spin-up and spin- 
down electrons, and fx and fy are the non-equilibrium 
distributions of the transverse spin components. The spin 
accumulation is defined as /x = deTr[o-/] 10]. The 
spin accumulation in the nonmagnetic layer is defined in 
a similar way. 

The spin accumulation in a nonmagnetic layer, /xn, 
obeys the diffusion equation [ll|, and is expressed as a 
linear combination of exp(±2:/Asd(N))j where Ag^jN) is the 
spin diffusion length. The spin current in a nonmagnetic 
layer is given by 

where S is the cross section area of the system, ctn is 
the conductivity and e is the absolute value of the elec- 
tron charge. The spin current in the N3 layer is equal to 
_jN3^F2 X — L2 + d2 because of the continuity of the 
spin current, and vanishes at x = L2 + d2 + L3 (see Fig. 
[T]). Using above boundary conditions, we obtain the spin 
accumulation in the N3 layer. 

The longitudinal spin current in a ferromagnetic layer 
is given by 

(m-lf)m = -^|^(crVF-CTVF)ni ' (4) 



where '''^ 



dx 2e2 

dep'^^^ is the electro-chemical poten- 



tial for the spin-up (spin-down) electrons and a"^^^^ is 
the conductivity of spin-up (spin-down) electrons. The 
polarization of spin-dependent conductivity is defined as 
f3 = (ct^ — a^)/{a^ + (J^). The longitudinal spin current 
in the F2 layer is equal to m2 • 1^2 at x = ^2 and 
— m2 • aX X = L2 + d2 because of the continuity 

of the spin current. Solving the diffusion equation [ll| 
with the above boundary conditions, we obtain longitu- 
dinal spin accumulation in the F2 layer. The longitudinal 
spin accumulation is expressed as a linear combination of 
exp(±x/As(j(pj^)), where Asd(pL) is the longitudinal spin 
diffusion length. 

We assume that the transverse spin accumulation 
obeys the following equation 0]: 
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^/^F = ^A'F X m+ ^2^^ fij , (5) 
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where Aj = ^/{D^ + D^)h/{2J) and Asd(FT) is the 
transverse spin diffusion length. Here J represents 
the strength of the exchange field. The transverse 
spin accumulation is expressed as a linear combina- 
tion of exp(±a://+) and exp(±x/Z_), where l/l± — 

y^(-'-/'^sd(FT)) "F (i/'^j)- Therefore, we define the pene- 
tration depth of the transverse spin current At by 



— = Rc 
At 



(6) 
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The transverse spin current in a ferromagnetic layer is 
expressed as 



m X (I^ X m) 



dx 2e 



T 

2 f^F 



(7) 



where cr^^ = (l/2)(crV(l + + ^7(1 " /?'))• Here 
(3' = - D^)/{D^ + D^) is the polarization of the 
spin-dependent diffusion constants, and Q- For 
simiplicty, we assume that (3 — [3 . The transverse spin 
accumulation in a ferromagnetic layer is obtained by solv- 
ing Eq. ([5]) with boundary conditions satisfying the con- 
tinuity of the spin current at the N/F interface. 

Solving the diffusion equations of the spin accumula- 
tions of the N3 and F2 layers, the backflow at the N2 /F2 
interface is re-written as 
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(8) ■ 

where the conductance g*^^) is given in Ref. [lo| . and de- 
pends on the ratio \sd(Fi^)- Similarly, the renormal- 
ized mixing conductances, gl\(Y^)^ depend on the ratio 
d2/l+{F^)- If the thickness of the N3 layer is thin enough 
compared to its spin diffusion length, gWfY^) ^'^^ given 
by 



where A 
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K2^ and Ki and K2 are given 
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where gt/S — h/2e'^ pp^l^ and pp.-^ is the resistivity of 
the F2 layer. The mixing conductance of the Fi layer 
in Eqs. H]) and ([2]) is also replaced by the renormalized 
conductance. 

We assume that spin-flip scattering in the N2 layer 
is so weak that we can neglect the spatial variation 
of the spin current in the N2 layer. Then we have 
jpump _ jNi^Fi ^ jNs^Fa^ a,nd the spin accumulation 

in the N2 layer can be determined The torque act- 
ing on the magnetization of the Fi layer is given by 
mi X {(IP™P - lN2^Fi) X j^^i^ ^^^^^ yigyg ^Yie fol- 
lowing modified Landau-Lifshitz-Gilbert (LLG) equation 

Mm ■■ 

—J— = -7offmi X Beff H [ao + a )mi x — — , (12) 
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FIG. 2: The dependences of the line width of the FMR power 
absorption spectra, AB, on the thickness of the F2 layer, d2. 
Materials of the F2 layer are (a) Ni8oFe2o, (b) Co75Fe25 and 
(c) (Co5oF5o)8oB2o, respectively. The filled circles represent 
experimental data and the solid lines are fit to the experimen- 
tal data according to the theory with the finite penetration 
depth of the transverse spin current At. The dotted lines rep- 
resent the case of At = 0. 



where Bos is the effective magnetic field, 7 is the gy- 
romagnetic ratio, ao is the Gilbert damping constant 
intrinsic to the ferromagnetic metal, and a' is the en- 
hancement of the Gilbert damping constant due to spin 
pumping. The Gilbert damping constant is related to the 
line width of the FMR absorption spectrum via [l^ 



AB = ABo 



2lu 



(13) 



where lu = 27r/ is the angular velocity of the oscillating 
magnetic field. We notice that the effects of the Ni and 
N3 layers are quite small because, as mentioned below, 
the thickness of these layers are thin enough compared 
to its spin diffusion length in our experiments. Assuming 



that gl^ g}^ [H, in the limit of 6* ^ 0, we find 



AB - ABo 
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where 



r(F,)^ 



1,2) is the real part of the renormal- 



ized mixing conductance of the i-th ferromagnetic layer. 
We should note that if we neglect the transverse spin 
accumulation in the ferromagnetic layer the mixing con- 
ductances are not renormalized, and that the line width 
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Ai? does not depend on the thickness of the F2 layer [l^l ■ 
This is due to the fact that the dominant component of 
the pumped spin current is perpendicular to the magne- 
tization vector m2 of the F2 layer in our experiment. 

We performed FMR experiments on the three differ- 
ent Ni/Fi/ N2/F2/N3 five-layer systems shown in Fig. [T] 
1J| . Nonmagnetic layers are made of Cu. The combina- 
tions of the ferromagnetic layers (Fi,F2) of each system 
are (a) (CoFe,Py), (b) (Py,CoFe) and (c) (CoFe,CoFeB). 
The samples were deposited on Corning 1737 glass sub- 
strates using an rf magnetron sputtering system in an ul- 
trahigh vacuum below 4xlO~^ [Pa] and cut to 5 [mm^]. 
The Ar pressure during deposition was 0.077 [Pa]. The 
thickness of all Cu layers are 5 [nm]. The thickness of Fi 
layers is 5 [nm] for sample (a) and (b), and 10 [nm] for 
sample (c). The FMR measurements were carried out 
using an X-band microwave source (/ = 9.4[GHz]) at 
room temperature. The microwave power, modulation 
frequency, and modulation field are 1 [mW], 10 [kHz], 
and 0.1 [mT], respectively. The precession angles of all 
samples are estimated to be 1 [deg]. The resistivity pp 
of Py, CoFe and CoFeB are 241 [17nm], 94 [Onm] and 
1252 [f7nm], respectively. The magnetizations (AirM) of 
Py and CoFe are 0.76 [T] and 2.1 [T], respectively. The 
gyromagnetic ratio is 1.8467 x 10^^ [Hz/T] for all systems 

In Figs. [11(a), (b) and (c), the measured line widths of 
the FMR absorption spectra AB are plotted with filled 
circles against the thickness of the F2 layer, d2- The solid 
lines are fit to the experimental data according to the 
theory with the finite penetration depth of the transverse 
spin current At . The dotted lines represent the calculated 
AB in the case of At = [I^. 

Parameters other than At are determined as follows. 
The mixing conductances per unit area of the combina- 



tions {gl^p^^/S, gl^p^s^/S) are assumed to be (a) (48.0, 

38.0), (b) (15.2, 17.0) and (c) (48.0, 128.0) [nm-^]. Al- 
though these values are determined by fitting, they have 



good agreement with the ab initio caluclations P. For 
simplicity, we assume that 4^ = tV where values oitiyS 
of Py, CoFe and CoFeB are taken to be 4.0 [nrn'^], 6.0 
[nm~^] and 0.8 [nm~^], respectively. The longitudinal 
spin diffusion lengths are 5.5 [nm] for Py and 12 [nm] for 
CoFe and CoFeB, respectively [lH[3- The polarizations 
of conductance (3 are 0.73 for Py, 0.65 for CoFe and 0.56 
for CoFeB, respectively p5l[l^.Tl7j. The transverse spin 

diffusion lengths are given by Asd(FT) = Kd(Fi,)/\/^ ~ 1^^ 



We take qI^JS = 1.0 [nm'^]. 



2gTTgU/(5TT_ 



20.0 [nm~'^] [lO| for all systems; these are not impor- 
tant parameters for fitting the experimental results. The 
spin diffusion length and resistivity of Cu are taken to 
be 500 [nm] and 21 [Onm] [l^. The obtained values of 
At are 3.7 [nm] for Py, 2.5 [nm] for CoFe and 12.0 [nm] 
for CoFeB, respectively. Our results agree quite well with 
the prediction based on the Boltzmann theory of electron 
transport [3]. 

In conclusion, we analyzed spin pumping in 
Cu/CoFe/Cu/Py/Cu, Cu/Py/Cu/CoFe/Cu and 
Cu/CoFe/Cu/CoFeB/Cu five-layer systems both 
theoretically and experimentally. We showed that the 
enhancement of the Gilbert dumping constant due to 
spin pumping depends on the ratio of the penetration 
depth At and the thickness of the ferromagnetic layers, 
which is not in resonance with the oscillating magnetic 
field. We measured the line widths of FMR absorption 
spectra, which are closely related to the Gilbert dump- 
ing constant. Analyzing the experimental results, we 
showed that the penetration depths of the transverse 
spin current in Py, CoFe and CoFeB are 3.7 [nm], 2.5 
[nm] and 12.0 [nm], respectively. Our results support 
the Boltzmann theory of transverse spin current [7|. 
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